Critical scalings and jamming in thermal colloidal systems 
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During the jamming of thermal colloids, the first peak of the pair distribution function shows a 
maximum height gY^'^^. We find that gY^'^'^ is accompanied by significant change of material properties 
and thus signifies the transition from unjammed to jammed glasses. The scaling laws at gf^^ lead 
to scaling collapse of structural and thermodynamic quantities, indicating the criticality of the 
T = jamming transition. The physical significance of 3™^" is highlighted by its coincidence with 
the equality of the kinetic and potential energy and the maximum fluctuation of the coordination 
number. In jammed glasses, we find the strong coupling between the isostaticity and flattening of 
the density of vibrational states at the isostatic temperature scaled well with the compression. 
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As described by the seminal jamming phase diagram 
[H , a packing of frictionless spheres interacting via repul- 
sion undergoes the jamming transition at a critical point 
denoted as 'J', which generalizes the jamming of col- 
loids and granular materials [3-[|. At r = 0, marginally 
jammed solids near Point J show particular critical scal- 
ings and length scales [6l-[TT| . Although recent theoretical 
work suggests that the jamming transition does not di- 
rectly link to the glass transition in the hard sphere limit 
[3. [T2| - [l3 | ■ the T = jamming transition still sheds some 
light on understanding the nature of glasses and amor- 
phous solids [ill, mi • 

Like the long-standing glass transition problem [l6l . 
[iTj , the structural similarities and absence of long range 
order during the jamming transition obscure proper order 
parameters to show diverging static correlation length, so 
the nature of the jamming transition still remains elusive. 
Only recently, a structural vestige of the T = jamming 
transition has been observed in experiments and simu- 
lations of thermal colloidal systems 0, [l^- At Point J, 
particles are just in contact, which induces a 5— function 
of the first peak of the pair distribution function g(r) 
[UliOl. At r > 0, 31, the height of the first peak of g{r), 
reaches a maximum value g^^^^ at a temperature depen- 
dent volume fraction (pressure). In the T = limit, g™^^ 
diverges at Point J 

It remains an open question whether the emergence of 
^max jg merely the thermal vestige of the T = jamming 
transition. If not, does it imply any physical significance 
in the formation of amorphous solids? What features of 
the r = jamming transition would persist at T > 
and in what manner? In this letter, we provide strong 
evidence via molecular dynamics simulations at low tem- 
peratures to confirm that the emergence of 5™^'^ contains 
rich and important physics. It marks the transition from 
unjammed to jammed glasses by showing critical scalings, 
recovering of some T = jamming features, and remark- 
able signs in the fiuctuations and energy competition. 
In jammed glasses, we find a novel isostatic temperature 
Tj at which the isostaticity and flattening of the density 



of vibrational states coincide. The isostatic temperature 
is scaled well with the compression. Interestingly, the 
onset frequency of the plateau in the density of vibra- 
tional states at Ti shows similar pressure dependence to 
marginally jammed solids at T = 0. 

Our systems are three-dimensional cubic boxes con- 
sisting of = 1000 spheres with the same mass m 0. 
Half of the spheres have a diameter of a, while the other 
half have a diameter of lAa, to avoid crystallization. 
Periodic boundary conditions are applied in all direc- 
tions. Particles i and j interact via repulsive potential 
Vij = ^ (1 — i^ij/'^ij)" when their separation is less 
than the sum of their radii crij , and zero otherwise. In this 
letter, we only show the results for harmonic repulsion 
(a — 2). The results for Hertzian repulsion (a = 5/2) 
are shown in the Supplementary Information. We set the 
units of mass, energy, and length to be m, e, and ct, and 
the Boltzmann constant fc^ = 1 . The frequency is in the 
units of ^fTfma^ . We perform molecular dynamics simu- 
lations at constant temperature and pressure using Gear 
predictor-corrector algorithm. The density of vibrational 
states D(uj) is obtained from the Fourier transform of the 
velocity correlation function C(t) ~ {v{t) ■ v{0)) (2l| : 
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where (.) denotes the time and particle average. 

Figure [T] provides sufficient information to point out 
the physical significance of the structural signature gf^'^^. 
The left panels of Fig. [1] indicate that along with the 
emergence of g™^^ at a crossover pressure pj multiple 
quantities undergo qualitative changes. In a wide range 
of temperatures, both the volume fraction difference from 
Point J (f) — (j)c and potential energy per particle V show 
distinct pressure dependence on two sides of pj . At p > 
Pj, some of the critical scalings in marginally jammed 
solids at T = are recovered: 
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FIG. 1: (color online) Left panels: Pressure dependence of the 
height of the first peak of the pair distribution function gi, 
volume fraction difference from the T = jamming transition 
(f> — <f)c, potential energy per particle V, and mean square 
fluctuation of the coordination number Sz'^ for systems with 
harmonic {a = 2) repulsion. The black solid lines show the 
scalings in marginally jammed solids at T = 0: (j) ^ i'c ^ 
pI<\? and V ~ (p/(j!>^)^. The dashed lines mark the locations 
of f/f^''. Right panels: Scaling collapse of the quantities in 
the left panels. The data are measured at T = 10"* (red 
circles), 10~^ (blue squares), 4 x 10^^ (maroon pluses), 10~® 
(green triangles), 4 x 10~^ (magenta stars), and 10^^ (orange 
diamonds) . 



where the critical volume fraction at Point J c/)c = 0.649 ± 
0.002 for both harmonic (a = 2) and Hertzian (a — 2.5) 
repulsions. We correct the widely accepted scaling p ^ 
((/) — (/)c)"~^ to Eq. ([2]) which can be simply derived 
from Eq. ^ and -p = 0^^, because away from 0^ the 
variation of must be considered to interpret the data 
correctly. At p < pj, however, the jamming-like scalings 
break down, indicating the notable change of the material 
properties across . 

In addition to Eqs. ([1]) and ([3]), the average coordina- 
tion number per particle z is scaled well with the com- 
pression at T^O: z-z^^ (0-0c)^^^ ~ (p/(/«2)i/2(a-i)^ 
where Zc = 2d is the isostatic value with d the dimension 
of space. Isostaticity (z = Zc) is one of the most special 
features of the T = jamming transition at Point J. It 
controls the unusual vibrational properties of marginally 
jammed sohds and realistic glasses [l-H, [lH. Its role in 
vibrations of jammed glasses at T > will be discussed 
later. At T > 0, the thermal motion breaks particle 
contacts frequently, so the T = scaling of z no longer 
exists in the vicinity of g^^^^. We observe that Zj, the 
coordination number at pj is less than Zc and gradually 



decreases with decreasing the temperature. Meanwhile, 
the slope {^)zj grows rapidly to infinity in the T = 
limit. This is consistent with the discontinuous jump of 
the coordination number from to Zc at Point J @. 

Although the coordination number does not maintain 
the T = scaling at p > pj, its fluctuation Sz"^ is inter- 
estingly correlated to g™''^ by peaking at Pj, as shown 
in the left bottom panel of Fig. [1] It is natural to imag- 
ine that the fluctuation of particle contact is smaller in 
more solid-like systems. The peak in 5z'^ indicates that 
systems at p > pj is more rigid than those at p < pj. 
Furthermore, the peak value of Sz"^ is almost indepen- 
dent on the temperature at low temperatures, implying 
the singularity of Point J, since Sz"^ = at T = 0. 

A close look at the potential energy unveils another 
interesting phenomenon: V k, dX pj. At p > pj, the 
potential energy dominates the kinetic energy, so that 
the thermal motion is not strong enough to cause sig- 
nificant configuration change of jammed states and the 
T = scalings thus hold. At p < pj, the kinetic energy 
wins. The thermal motion has noticeable effects on the 
system performance. Therefore, the emergence of g™^^ 
and the corresponding change of material properties do 
not happen by chance, but have explicit physical origins. 

At T = 0, the height of the first peak of the pair dis- 
tribution function gi is inversely scaled with cj) — (pc when 
the system is compressed away from Point J, because 
gi is inversely proportional to the particle overlap which 
grows linearly with (j) — (j)c @j [3 ■ It has been shown that 
Tj , the temperature at pj is scaled well with (j) — (j)^ 0] : 

T,^(^-0.)"^(p/0^)"/^"-^^ (4) 

Neglecting the thermal effects on gi from T = to , we 
can easily obtain g™^^ ^ (0 — (j)c)~^ ~ j^-i/a^ However, 
thermal and nonlinear effects still cause visible deviation 
from the prediction. As shown in the Supplementary 
Information, 

gr''-iP-<i^cr'<-{p/^T"^^"-'\ (5) 

where 7 = 0.90 ± 0.02. 

Equations ((H) — (H]) imply the criticality of Point J 
viewed in thermal colloidal systems. By scaling all the 
structural and thermodynamic quantities in the left pan- 
els of Fig. [1] using Eqs. (H]) — (O, we obtain excellent scal- 
ing collapse of the original low-temperature data. As 
shown in the right panels of Fig. [U there are continuous 
scaling functions covering data on both sides of pj : 

where ^ denotes gi, (j> — 4>c, V , and Sz"^ with ~ —7/0;, 
l/a, 1, and respectively. During the writing of this 
letter, we notice that very recent approaches based on 
mean field theory has predicted similar scaling collapse 
[2^ . [2^ . Here we explicitly obtain Eq. ^ directly from 



3 




"(b) 










FIG. 2: (color online) (a) Density of vibrational states -D(i^) 
for systems with harmonic (a = 2) repulsion at T = 10"'''. 
The pressures are 0.0008, 0.001, 0.002, 0.003, 0.004, 0.005, 
and 0.009 from the black to violet curves. The crossover pres- 
sure pj is around 0.0015. (b) Slope of D{uj) in the interme- 
diate frequency regime versus the coordination number z for 
T = 10~* (red circles), 10"^ (blue squares), 10"^ ( green tri- 
angles), and 10"^ (orange diamonds). The dashed line shows 
where the flattening of D{ij) happens, (c) Frequency uj* close 
to the onset of the plateau in D{ljj) with D{u)*) ~ 0.66 mea- 
sured at 2 = 2c. The dashed line shows the T = Q scaling. 



the scalings present at pj without making any assump- 
tions. Furthermore, the scahng collapse of 5z^ is not able 
to predict from the mean field theory. Our study here 
gives exphcit physical origin and meaning of the scaling 
collapse, and highlights further the physical significance 
of g^^^. Since the scaling collapse is simply derived from 
the scahng laws at pj which approaches Point J in the 
T = and p = Q limit, our finding here strongly supports 
the criticality of Point J. The scaling collapse conveys im- 
portant information of the low-temperature properties of 
the model systems in the vicinity of Point J, not just 
trivially refiecting the behaviors at T = 0. 

From the above discussions about Fig. [U we can con- 
clude that the emergence of g^'^^ is not trivially the ther- 
mal vestige of the T = jamming transition. It is plau- 
sible to state that g'f'^^ is the structural signature of the 
jamming transition from unjammed to jammed glasses 
at T > 0. 

As mentioned above, the T = scaling of the coordi- 
nation number z is violated in jammed glasses near pj. 
At Pj, z < Zc although <f> > (pc- At fixed temperature, 
z increases with the pressure, so isostaticity (z = z^) 
would occur at p > pj. At T = 0, isostaticity determines 
the special vibrational properties of marginally jammed 
solids 0, BEll- For instance, the density of vibrational 
states D{uj) possesses a low-frequency plateau at Point 
J, which shrinks with increasing (f) — (f>c and eventually 
disappears 0]. The characteristic frequency of the onset 
of the plateau is linearly scaled with z — Zc for jammed 
systems with harmonic repulsion 0, Q • It is thus inter- 
esting to know if there arc any connections between the 



isostaticity and vibration in jammed glasses at T > 0. 

Figurema) shows the evolution of D{uj) measured from 
Eq. (dl with the change of the pressure at fixed temper- 
ature. Notice that all the systems shown here are glassy 
with the relaxation time much longer than the simulation 
time window. At low pressures, there is a low-frequency 
peak which grows up with decreasing the pressure. The 
same trend has also been observed in hard sphere systems 
[2^ . This low- frequency peak indicates the aggregation 
of soft modes, so the system with a higher peak would be 
less stable. At pj, this soft mode peak is still pronounced, 
in consistent with the fact that z < Zc- At p > pj, the 
peak is progressively suppressed with the pressure. Inter- 
estingly, a low-frequency plateau is eventually formed in 
D{uj), which is also shown in a recent experiment of ther- 
mal colloids [2^. Afterwards, the pressure dependence of 
D{Ld) looks similar to that at T = 0. 

We estimate the average slope of D{u!) in the inter- 
mediate frequency regime at different temperatures, as 
shown in Fig. ^h) . It is interesting that the flattening of 
D{iL!) happens approximately when z = Zc at the isostatic 
temperature Ti. This striking finding indicates that the 
flattening of D{Ld) is strongly coupled to the isostaticity, 
not only at Point J. but also in thermal colloidal systems 
with repulsion. At Ti, although z = Zc and there is a 
plateau in D{uj), the onset frequency uj* of the plateau 
still relies on the pressure. Figure [5Jc) shows that the 
T = scaling of jammed solids with harmonic repulsion, 
- (0 - (jicY/"^ - {pl<i?fl'^ 0,il, is stih valid at Ti. 

Figure [5] unveils the complicated and puzzling effects 
of the thermal motion on the vibration. Notice that the 
reference state at T^, i.e. the configuration with time- 
averaged particle positions, still has a coordination num- 
ber larger than Zc, although in time average z = Zc. Iso- 
staticity induced by thermal fluctuations only controls 
the flattening of T){bS). The compression, or equivalently 
the reference state, may still determine some vibrational 
properties of jammed glasses. Here we have the first at- 
tack on this interesting issue. Further studies are needed 
to better understand the thermal effects. 

The isostatic temperature Ti is scaled well with the 
compression: 



T 



a+1 
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(7) 



This scaling may be understandable from the competi- 
tion between kinetic energy and potential energy to break 
extra contacts of the reference states beyond isostaticity. 
Together with Eq. ^ and a recent observation of the lin- 
ear pressure dependence of the glass transition temper- 
ature at which the relaxation time or viscosity diverges, 
Tg ^ p [26j . we are able to plot the phase diagram in 
Fig. El to clarify the connections between different transi- 
tions. At low temperatures, the glass transition, jamming 
transition, and isostaticity happen in sequence when the 
pressure (volume fraction) increases. We have already 
known that there is a gap between the glass transition 
and jamming transition [3, [27| . Here we introduce the 
new isostaticity line which converges to Point J with the 
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FIG. 3: (color online) Phase diagram of systems with har- 
monic (q — 2) repulsion. The glass transition temperature 
(blue diamonds), jamming transition temperature (black cir- 
cles), and isostatic temperature (red squares) are plotted to- 
gether as a function of p/4>'^ ■ The lines are power-law fits 
to the data: Tg ~ (blue dashed), Tj ~ {p/(t>^f (black 

solid), and Ti ~ {j>/(f>^)^ (red dot-dashed). 

jamming transition in the T = limit, but likely on the 
two limits of Point J since z is different on these two lines 
and discontinuous at Point J. Figure |3] demonstrates the 



complexity in the formation of amorphous solids. Even in 
glasses, there are complicated thermodynamic and vibra- 
tional behaviors beyond our expectation, which requires 
intensive studies to sort through. 

The findings discussed in this letter are limited to low 
temperatures and pressures where the influence of Point 
J is still sensible. As shown in Fig. [3l if all the critical 
scalings hold at high temperatures, the three lines even- 
tually get across. Recent studies have shown that the 
glass transition temperature drops with increasing the 
pressure at sufficiently high pressures [28| . so nothing at 
high temperatures and pressures is predictable from the 
low-temperature and low-pressure properties. Beyond 
the regime concerned here, we will enter a new territory 
with completely new deep jamming pictures (29j . 
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FIG. 4: Temperature dependence of the maximum value of 
the first pealt of tlie pair distribution function, g™^". The red 
circles and blue squares are for systems with harmonic (a = 2) 
and Hertzian {a = 2.5) repulsions, respectively. The red and 
blue dashed lines show the power law scalings: gf ~ y-o.46 
and 35"'"' ~ T'^-'^^ (g^'"' ~ respectively. 




FIG. 5: Scaling collapse of the height of the first peak of the 
pair distribution function gi, volume fraction difference from 
the r = jamming transition (f) — (jic, potential energy per 
particle V, and mean square fluctuation of the coordination 
number 5z^ for systems with Hertzian (a = 2.5) repulsion. 




p/(^' 

FIG. 6; Phase diagram of systems with Hertzian {a = 2.5) 
repulsion. The glass transition temperature (blue diamonds), 
jamming transition temperature (black circles), and isostatic 
temperature (red squares) are plotted together as a func- 
tion of p/4>^. The lines are power- law fits to the data: 
Tg ~ p/cfp (blue dashed), Tj ~ {p/'p^f^''^ (black sohd), and 
- {j)/(l>'^y^''^ (red dot-dashed). 
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